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Abstract. This pedagogical review aims at presenting the fundamental aspects of 
the theory of inflationary cosmological perturbations of quantum-mechanical origin. 
The analogy with the well-known Schwinger effect is discussed in detail and a sys- 
tematic comparison of the two physical phenomena is carried out. In particular, it 
is demonstrated that the two underlying formalisms differ only up to an irrelevant 
canonical transformation. Hence, the basic physical mechanisms at play are simi- 
lar in both cases and can be reduced to the quantization of a parametric oscillator 
leading to particle creation due to the interaction with a classical source: pair pro- 
duction in vacuum is therefore equivalent to the appearance of a growing mode for 
the cosmological fluctuations. The only difference lies in the nature of the source: 
an electric field in the case of the Schwinger effect and the gravitational field in 
the case of inflationary perturbations. Although, in the laboratory, it is notoriously 
difficult to produce an electric field such that pairs extracted from the vacuum can 
be detected, the gravitational field in the early universe can be strong enough to 
lead to observable effects that ultimately reveal themselves as temperature fluctua- 
tions in the Cosmic Microwave Background. Finally, the question of how quantum 
cosmological perturbations can be considered as classical is discussed at the end of 
the article. 



1 Introduction 

The scenario of inflation was invented in order to solve puzzling issues associ- 
ated with the standard hot Big Bang theory [II H]. Soon after its advent, it was 
realized that this scenario also contains a remarkable extra bonus: it gives a 
well-motivated mechanism for structure formation that leads to a nearly scale- 
invariant power spectrum [3,4, namely exactly what is needed in order to 
account for various astrophysical observations in a satisfactory way [5] . How- 
ever, this is not the only aspect that deserves to be stressed. Indeed, even 
from a fundamental point of view, this mechanism appears quite remarkable 
in the sense that it combines general relativity with quantum mechanics. The 
main purpose of this review article is to thoroughly discuss this aspect of the 
theory of inflationary cosmological perturbations. 
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This theory is in fact remarkable at two levels. Firstly, because it relies on 
the phenomenon of particle creation which is a non-trivial effect in quantum 
field theory. In this sense, it is equivalent to the well-known Schwinger [6] effect 
and this analogy will be made explicit in this paper. The basic ingredient is a 
quantum scalar field <P (in practice this is rather a fermionic field \P but, for 
simplicity, we will restrict ourselves to the case of a scalar field) interacting 
with a classical source, in the case of the Schwinger effect, an electric field 
E. The Schwinger effect has not yet been observed in the laboratory as it is 
difficult to produce an electric field with the required strength but there are 
prospects to do so, in particular at DESY with a Free Electron Laser (FEL) 
in the X-ray band [3 [9] but also at SLAC with the Linac Coherent Light 
Source (LCLS) [10] . Even if there is absolutely no reason to doubt the reality of 
the Schwinger effect, observing pair creation in the laboratory would clearly 
be a breakthrough and, in some sense, a verification of the corresponding 
inflationary mechanism. 

Secondly, the theory of cosmological perturbations is also remarkable for 
the following reason. In cosmology, what plays the role of the constant elec- 
tric field E [originating from a time-dependent potential vector A^{i)] is 
the background gravitational field, i.e. the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) scale factor a(t), and what plays the role of the quantum 
fermionic field <P(t, x) is the quantum perturbed metric 6g/j, v (t,x), that is to 
say the small inhomogeneous fluctuations of the gravitational field itself [IT] . 
In the early Universe, the gravitational field is quite strong, i.e. for instance 
H/m pl ~ 10~ 5 , where H is the Hubble parameter, and this is why the cos- 
mological version of the Schwinger effect can be efficient. From the previous 
considerations, it is also clear that, in some sense, the inflationary mecha- 
nism relies on quantum gravity which adds another interesting aspect to the 
problem. Of course, we only deal with linearized quantum gravity and this 
is why we do not have to face tricky questions associated with finiteness of 
quantum gravity and/or renormalization. More precisely, in the case of scalar 
perturbations, Sg^it, x) is replaced by the Mukhanov-Sasaki variable v(t, x) 
which is a combination of the Bardeen potential (the generalization of the 
Newtonian potential in general relativity) and of the small fluctuations in the 
inflaton field. For gravitational waves, the relevant quantity is hij(t,x), the 
transverse and traceless part of the perturbed metric. 

Let us notice that the two above-mentioned aspects are features of the 
theory of cosmological perturbations in general. The inflationary aspect is 
in fact not necessary in order to have particles creation: only a dynamical 
background is required. However, a quasi-exponential expansion is mandatory 
if one wants to obtain a power spectrum which is close to scale invariance as 
indicated by astrophysical observations. 

The fact that the inflationary mechanism for structure formation relies 
on general relativity and quantum mechanics also raises fundamental inter- 
pretational questions. In particular, the question of how classicality emerges 
is of special relevance in this context p~21 [13] . Indeed, the perturbations are 



The Cosmological Schwinger Effect 



3 



of quantum-mechanical origin but no astrophysical observations suggest any 
typically quantum-mechanical signature. Therefore, it is necessary to under- 
stand how the perturbations have become classical (and in which sense). This 
leads to very deep issues. For instance, if one invokes a mechanism based on 
the phenomenon of decoherence [14] , then one has to discuss what plays the 
role of the environment. This question is clearly non-trivial in the cosmologi- 
cal context. At the end of this review article, we will address these questions 
using the Wigner function [15] as a tool to understand when a system can be 
considered as classical. 

This paper can be viewed as the third of a series on the inflationary the- 
ory, the two first ones being Refs. [TBI HZ] . The topics developed in those 
last two references will be supposed to be known and we will often refer to 
them. The present paper is organized as follows. In Sec. [2) we briefly review 
the Schwinger effect. In particular, we derive the rate of pair production in 
the Schrodinger functional approach and stress the importance of the Wentzel- 
Kramers-Brillouin (WKB) approximation as a method to choose a well-defined 
initial state. In Sec. El we quantize a free scalar field in a FLRW Universe and 
show that the basic physical phenomenon at play is equivalent to that respon- 
sible for the Schwinger effect, namely particle creation under the influence of 
a classical source. In particular, we demonstrate that, up to a canonical trans- 
formation, the underlying formalisms are the same. Roughly speaking, in both 
cases, one has to deal with parametric oscillators. The only difference between 
the two systems lies in the time dependence of the corresponding effective fre- 
quencies. In Sec. HI we argue that the equations obeyed by the cosmological 
perturbations (in particular during inflation) are equivalent to the equations 
of motion of a free scalar field. We emphasize that the relevant observable 
quantity is the two-point correlation function since it is directly linked to the 
Cosmic Microwave Background (CMB) temperature fluctuations. Finally, as 
mentioned above, in Sec. El we address the question of the classicality of the 
cosmological perturbations. 

2 The Schwinger Effect 
2.1 General Formalism 

The action of a complex (charged) scalar field interacting with an electromag- 
netic field is given by 



where rf 1 ^ is the flat (Minkowski) space-time metric with signature (—,+,+,+) 
and where the covariant derivative can be expressed as 




(1) 



V a $ = d a <P + iqA a <P , 



(2) 
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q being the charge of the field. The quantity to represents the mass of the 
scalar particle. Assuming the following configuration for the vector potential 
= (0, 0, 0, —Et), where E is the magnitude of the static electric field aligned 
along the z direction (by convention), one obtains the following equation of 
motion 

dt$ - d i d l $ + 2iqEtd z <P + q 2 E 2 t 2 <P + m 2 $ = . (3) 

It turns out to be more convenient to Fourier transform the field since this 
allows us to study the evolution of the system mode by mode. For this purpose, 
one decomposes the field according to 

<i>(t,x) = -^Jd 3 k<l> k (t)e ik -*. (4) 

In the above expression, <P k {t) is the time-dependent Fourier amplitude of the 
mode characterized by the wave- vector fc. Inserting the Fourier transform (j4]) 
into Eq. (Q]), the action of the system takes the form 

S = - [ dt / dfc \-& k $i + (k 2 - 2qEk z t + q 2 E 2 t 2 + to 2 ) $ k $* k 1 , (5) 

J JR 3 L - 1 

where a dot denotes a derivative with respect to time. The variation of this 
Lagrangian with respect to <P* k and leads to 

a r a r 

— = - (k 2 - 2qEk z t + q 2 E 2 t 2 + to 2 ) <P k , = p k = <P k , (6) 

s $ k s®% 

where p k is the conjugate momentum of the Fourier component of the field 
and C denotes the Lagrangian density in Fourier space. Using the two above 
formula, the Euler-Lagrange equation of motion reads 

<P k +iu 2 (k,t)$ k -0, (7) 

where the time dependent frequency uj(k, i) can be expressed as 

UJ 2 (k,t) = k 2 ~2qEk z t + q 2 E 2 t 2 + m 2 . (8) 

Eq. ([7]) is of course similar to the one one would have obtained by directly 
substituting Eq. (01 into Eq. © . It is the equation of motion of a parametric 
oscillator. Let us recall that a parametric oscillator is an harmonic oscillator 
whose frequency depends on time. A typical example is a pendulum with a 
varying length. 

Let us now pass to the Hamiltonian formalism. The Hamiltonian is ob- 
tained from the Lagrangian by a standard Legendre transformation and can 
be expressed as 



H = f dk ( Pk S* k + P * k & k -£)= f dfc [p k p* k + uj 2 (fc, t) $ k $* k ] 

Jw 3 V ' JR 3 



(9) 
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For the following considerations, it turns out to be convenient to also work 
with real variables instead of the complex Therefore, we now introduce 
the definitions 



Then, the Hamiltonian can be written as 
H= [ dfe 



~ (Pi) 2 + \J (k,t) {nf + \ (Pif + \J (M) ($l) 2 



(10) 



(11) 



One recognizes the Hamiltonian of a collection of parametric oscillators with 
a time-dependent frequency given by Eq. ((8|). Again, one can check that the 
Hamilton equations deduced from Eqs. ([9]) and (fTTjl lead to an equation of 
motion similar to the one already derived before, namely Eq. ([7]). 



2.2 Quantization 



Our next step is to describe the quantization of the system. More precisely, 
the complex scalar field is quantized while the gauge field remains classical. 
Therefore, we have to deal with the interaction of a quantum field with a clas- 
sical source. Quantization is achieved by requiring the following commutation 
relations (a hat symbol is put on letters denoting operators) 



(k-p) , 



iS^ (k - p) 



(12) 



We choose to work in the Schrodinger picture where the states are time- 
dependent and the operators constant. The above commutation relations ad- 
mit the following representation 



(13) 



The state of the system is described by a functional of the scalar field, 
fy[<P(t : x)] (in the present context, ^ is the field functional and has clearly 
nothing to do with the fermionic field mentioned before), which can also be 
viewed as a function of an infinite number of variables, namely the values of 
<P at each point in space. Alternatively, one can also consider this functional 
as a function of the infinite number of Fourier components of the field and 
write 



(14) 



In the above equation, n represents the number of modes that, in the interme- 
diate calculations, it is useful to keep finite (for instance, if we imagine that 
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the field lives in a finite box). However, at the end, we will always consider 
the continuous case and take the limit n — > +00. 

In the framework described before, the Schrodinger equation is a functional 
differential equation. However, the Hamiltonian takes the form of an infinite 
sum over fc, see for instance Eqs. ([9]) and (fTTj) . without any interaction between 
different modes. As a consequence, each mode evolves independently and the 
corresponding Hamiltonian is represented by an ordinary differential operator. 
Explicitly, one has 

1 <9 2 <f 1 o 9 1 d 2 >P 

H k * = (flj + H k ) V = + ^ 2 (fc, t) {$lf* - 



2 8(4%)' 2 2 d (& h y 



1 



y(k,t)(& k ) 2 <P, (15) 

where the frequency uj has been given by Eq. ([5]) . 

Let us now consider the ground state of the system described before. We 
will discuss the choice of the initial conditions and the meaning of the vacuum 
state in the following but, as is well-known, it is given by a Gaussian state 

(t,$l) = N k (t)e- n "^T , * l k {t,$ k ) =7V fc (;)e- fifc(t) K) 2 , (16) 

where N k (t) and f2 k (t) are functions that can be determined using the 
Schrodinger equation id^f = H^. This leads to 

& = n k , h k = -unl + 1 -lu 2 (k, t) . (17) 

I\ k z 

The equation for the complex quantity f2 k is a non-linear Ricatti equation. 
When a particular solution is known, the general solution can be found by 
means of two successive quadratures [TS] . But this non-linear first order differ- 
ential equation can also be transformed into a linear second order differential 
equation. It turns out that this last one is exactly the equation for the Fourier 
mode function, Eq. ([7]). Therefore, the solutions to Eqs. (| 1 T|) read 

/23W4\ 1/4 ^ ifh 



n-y-ir) ■ °> = - 2 t- < 18 » 

where f k obeys the equation f k + oj 2 (k,t)f k — 0, that is to say, as already 
mentioned, the same equation as the Fourier component of the field, namely 
Eq. ([7])- The quantity N k is obtained by normalizing the wave-function. One 
can check that this leads to an equation consistent with the first formula in 
Eqs. (fl7|) . Therefore, one obtains that the ground quantum state of the field 
is given by 



k 



^23W4y /2 c -, 2fc(f) [ K) VK) 2 



(19) 
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Once a particular solution for the mode function has been singled out, the 
functions TV*, and flf~ , and hence the wave- function of the field, are completely 
specified. 

One can now use the above-mentioned state in order to calculate the am- 
plitude associated with the transition between two states <Zq and \&2- It is 
defined by 



0W = /nd*8<^(^r* 



1/2 fmo^k 



1/2 



xe 



EpK, P +^2, P )K) 2 



(20) 



from which, after having performed the Gaussian integration, one deduces 
that 



|<!Pi|!Pa 



= dct 



^2,k ) ( £l\,h + 



2.fc 



(21) 



At this point, one has to use the specific form of f2 k , in particular its expression 
given by Eq. (JTBJ) in terms of the function . One obtains 



det 



fl,kh,k ~ f2,kf* tk ) ( /2,fe/l,fe ~ /l,fe/2,fc 



(22) 



In this formula /i^ is the mode function for the initial state while / 2i fc is the 
same quantity but for the final state. As usual, one can always expand /2,fc 
over the basis (/i,fe, /i*.) and write 



h,k = Q!fc/i,fc + Pkf* k 



(23) 



Then, using the fact that the Wronskian fi,kf* k — fikfi,k is a conserved 
quantity (as can be easily checked by differentiating it and using the equation 
satisfied by f k ), one arrives at [It?] 



= dct 



1 



oik 



(24) 



Therefore, the determination of the transition amplitude amounts to inte- 
grating the equation controlling the evolution of the mode function. Once 
this is done, the coefficient a k is known and the quantity |(^i|!?2)| can be 
determined. In the next section, we discuss an explicit example. 



2.3 Particle creation 



We now use the formalism developed above in order to study the creation of 
quantum scalar particles due to interaction with a classical source. This is the 
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well-known Schwinger effect [8]. In the following, we will demonstrate that 
the inflationary mechanism for cosmological perturbations is exactly similar 
to the one discussed here, see also Refs. [213] . 

In order to determine the coefficient a k , one can proceed as follows. Let 
us use the dimensionless variable r = \JqEt — k z /^/qE. Then, the equation of 
motion for the Fourier component of the field, see Eq. (J7]), takes the form 

^ + (T + r 2 )^ = 0, (25) 

with T = {k\ + m 2 )/(qE). The quantity k ± is defined by k\ = k 2 — k 2 z — 
k 2 + k 2 (let us recall that we have chosen an electrical field aligned along the 
z-direction). Eq. (j25|) can be integrated exactly, see Eq. (9.255.2) of Ref. [21] 
and the solution can be expressed as 

$ k (r) = A k D_ (1+ir) /2 [(1 + i)r}+ B k D_ {1+ir)/2 [- (1 + i) r] , (26) 

where A k and B k are two constants fixed by the initial conditions and D p (z) 
is a parabolic cylinder function of order p. 

Despite the previous change of variable, Eq. (f2"5"|) has retained the form of 
an equation for a parametric oscillator but the frequency is now given by 



w(t) = Vr + r 2 . (27) 

This equation is well-suited to the WKB approximation. This approximation 
is not only useful to get an approximate form of the solution but can also be 
used in order to choose initial conditions that are well-motivated. Here, since 
we do already know the exact solution, it is clearly this last application we 
shall be concerned with. 

By definition, the WKB mode function (2uj)~ 1 / 2 e ±t $ " dT obeys the follow- 
ing equation of motion <P k + (co 2 — Q) <S> k = 0, where the quantity Q is defined 

by 

3 1 / dtuV 1 1 d 2 w 



Therefore, one sees that the WKB mode function is a good approximation to 
the actual one as soon as \Q/w 2 \ <C 1. This last condition defines the regime 
where the WKB approximation is valid. Let us compute this quantity in the 
case of Eq. (|25[) . Straightforward calculations lead to the following expression 



Q 



l l 



r 2 2(i + T 2 /ry 



5r 2 /r i 



2(i + r 2 /r) 



(29) 



The quantity \Q/uj 2 \ is represented in Fig. [TJ It is clear that, in the limits 
r/y/T — * ±oo, the WKB approximation is valid. This means that there exists 
a well-defined vacuum state (or adiabatic vacuum) in the "in" region, |0~), 
and in the "out" region, |0 + ). 



The Cosmological Schwinger Effect 



9 




Fig. 1. Evolution of the quantity \Q/lu 2 \ with time r in the case of the Schwinger 
effect. In the limit t/\/T — > ±00, \Q/ui 2 \ vanishes and the notion of adiabatic 
vacuum is available. 



When the WKB approximation is satisfied, an approximate solution of the 
mode functions is available and, as already briefly mentioned above, is given 

by 



$k (r) 



2uj 
a k <P 



exp 



wkb,fc 



;(0)d0 



(r) + /3 fc ^; kbjfe (r) 



exp 



+1 



uj{9)d9 



(30) 



where r; n i < is some arbitrary initial time. 

One can now use the WKB approximation in order to choose an initial 
state in the following way. We require that the system is in the adiabatic 
vacuum in the "in" region, |0 _ ) , that is to say when t/\^Y — > —00. Technically, 



this means that one has lim„ 



^wkb,fc or a k = 1 and (3 k 







(hence satisfying |a<k| — \[3 k \ =1)- This criterion completely specifies the 
coefficients A k and B k in Eq. and, as a consequence, also completely 
determines the coefficients a k and k m the "out" region (when r/\/T — > +00) 
that are needed in order to compute the transition amplitude. Around r ~ 0, 
see Fig. [IJ the WKB approximation is violated and we have particle creation. 
In the "out" region the vacuum is defined by |0 + ) and, therefore, the number 
of particles present in this region is measured by evaluating the amplitude 

(o-|o+). 

We now briefly explain how this can be done at the technical level. The 
phase can be computed exactly and reads 



10 



Jerome Martin 



In 



r 



T 



mi 

T 



(31) 



The arguments of the logarithms are always positive even if r is negative 
(hence the corresponding expression with ri n ; < is also meaningful). In the 
limit |t| T~ 1//2 ^> 1, the phase goes to 



jT w( 0)d^i(r|r| + |r ini | 2 ) + 



Y /It 



In 1-7-1 + In It-u 



and, therefore, the WKB mode function takes the form 



<Z> wkb , fc = -L e -^M+^)/3 | T |-i/a-*Hr/(ar) r , 



Then, in the limit r — > — oo, the exact solution given by Eq. 
expressed as 



(32) 

(33) 
can be 



2tt 



r(l-m/4 (1 _ ir (l-T)/2 ^2 e i| nnl | 2 /2 | T . n .|ir/2 $ 



wkb.fc 



Ake in(i + iT) /2 + Bfc ] (1 + ir ( l+ ir)/2 y2e-l-l 2 /2 ^i-ir/2 ^ kb ^ . 

(34) 



Since, as explained above, we choose the initial state such that a k = 1 and 
0k = 0, this amounts to requiring 



A k = 



r (2 + f) e ^(i-»r)/4 (1 _ ^(1-^/2 1 e - 4 |. lnl | 2 /2 | r . n .| 
v27r v 2 



iT/2 



Bfc = ^ fc e"( 1+ ^/ 2 . 



,(35) 
(36) 



Finally, one considers the behavior of the mode function in the limit r — > +oo 
and, using again the WKB mode function is this regime, one can find the 
corresponding coefficients a k and /3&. One obtains 



\a k \ 2 = l+e-* r , \p k Y 



(37) 



These expressions still satisfy \ctk\ — \pk\ — 1 as required. 

We have now reached our final goal and can return to the calculation of 
the determinant in Eq. ([24]) . Using the coefficient ak obtained above, one has 



|(o-|o+)| 



det 



1 



1 + e 



--rrT 



exp [— Trln (l 



wT 



)] 



(38) 



The evaluation of the trace is standard and leads to the well-known result 
first obtained by Schwinger in the early fifties [3] 
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VT 



OO 




(o-|o+> 



exp 



E 



{qE) 2 e 



nizm j (qE) 



(39) 



71=1 



The physical interpretation of this formula is clear. The argument of the ex- 
ponential gives the number of pairs created in the space-time volume VT 
due to the interaction of the quantum scalar field with the classical electric 
field. One can define a critical electric field (we have restored the fundamental 
constants) by 



which is such that the number of particles created is significant only if 
E 3> -Ecri- This condition can be understood by noting that this is just the re- 
quirement that the work performed by the force qE over the Compton length 
A = h/ (mc) is larger than the rest energy 2mc . In the case of pairs e + e~, the 
critical electric field is given by E crl ^ 1.3 x 10 18 V x m _1 . It is also interesting 
to remark that the dependence of |(0 _ |0 + )| in E is non-perturbative. This is 
one of the few example in quantum field theory where an exact result can be 
obtained (of course, this is not "full quantum theory" but rather "potential 
theory" since the radiative corrections to the Schwinger mechanism are not 
taken into account). 

We will see that the inflationary mechanism of production of cosmological 
perturbations is similar to the Schwinger mechanism. Therefore, observing this 
latter effect in the laboratory could be seen as an indication that we are on the 
right track as far as the inflationary mechanism is concerned. For instance, at 
DESY, there are plans to construct a Free Electron Laser (FEL) in the X-ray 
band which would effectively produce a very strong electric field and, hence, to 
observe the Schwinger mechanism [Tj [Sj ■ Unfortunately, even with a FEL, it is 
inconceivable to produce a static field with the required strength given present 
day technology. However, the situation is different for an oscillating electric 
field [22] (other configurations have been studied in Refs. [23] EH [25] ) and, in 
this case, it seems possible to extract pairs from the vacuum. This would also 
be a validation of the Schwinger mechanism since only the time-dependence 
of u>(k,t) is changed but not the other basic ingredients. It is also interesting 
to notice that, in the context of the inflationary theory, this case is in fact 
very similar to the reheating (26] [27] stage where the effective frequency of 
the perturbations is alternating due to the oscillations of the inflaton field at 
the bottom of its potential. 

To conclude this section, let us recall that the basic ingredient at play 
here is particle creation due to the interaction of a quantum field with a 
classical source. When the WKB approximation is valid, a well-defined notion 
of vacuum state exists, and when the WKB approximation is violated, particle 
creation occurs. We will see that, in the case of inflationary cosmological 
perturbations, exactly the same mechanism is available. 




qti 



(40) 
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3 Quantization of a Free Scalar Field in Curved 
Space-time 

Before considering inflation itself, let us now discuss the case of a free real 
scalar field in curved space-time since this is the simplest example which 
allows us to capture all the essential features of the theory of inflationary 
cosmological perturbations. 



3.1 General Formalism 



We consider the question of quantizing a (massless) scalar field in curved 
space-time. The starting point is the following action 

S=-J d^x^g^-d^d^ , (41) 

which, in a flat FLRW Universe whose metric is given by ds 2 = a 2 (r))(— dr] 2 + 
Sijdx l dx J ), rj being the conformal time, reads 

S=\J d 4 im 2 (?7) (<P' 2 - d'uK'N)/!') . (42) 

It follows immediately that the conjugate momentum to the scalar field can 
be expressed as 

n(r h x)=a 2 <P'(r ] ,x), (43) 

where a prime denotes a derivative with respect to conformal time. As before, 
it is convenient to Fourier expand the field ^(r/, x) over the basis of plane waves 
(therefore, we make explicit use of the fact that the space-like hyper-surfaces 
are flat). This gives 

^-iko^h 1 *^*- (44) 

We have chosen to re-scale the Fourier component with a factor 1/0(77) f° r 
future convenience. Since the scalar field is real, one has [i* k = fi-k- The next 
step consists in inserting the expression of <P(r], x) into the action (|42)l . This 
leads to 



S = - I dry / d^fc 



2 



a' 

Mfc Mfc + MfeMfc - 2— (Mfc^fe + Mfe Mfe) 

a 



a' 2 



a 2 ~ k ) (We/4 + MfcMfc) 



(45) 



Notice that the integral over the wave-numbers is calculated over half the 
space in order to sum over independent variables only |28j . This formula is 
similar to Eq. ([S]) for the case of the Schwinger effect (of course, in this last 



The Cosmological Schwinger Effect 



13 



case, we do not have <l> k = &*_ k since the field is charged and, hence, the 
integral is performed over all the momentum space). 

Equipped with the Lagrangian in the momentum space (which, in the 
following, as it was the case in the previous Section, we denote by C), one can 
check that it leads to the correct equation of motion. Since we have S£/Sfit = 
l/2[-2Hf/ k + 2(H 2 - k 2 )/J,' k ], the Euler-Lagrange equation d[SC/Sfi' h *]/dri - 
5£/6[i k — reproduces the correct equation of motion for the variable /Xfc, 
namely 

^±+^(k,r,)^ k = 0, (46) 

that is to say, again, the equation of a parametric oscillator, as in Eq. ||7J), but 
with a frequency now given by 



! (k,v) = k 2 --. (47) 



a 



This last formula should be compared with Eq. ([8|). In the case of the 
Schwinger effect, the frequency was time-dependent because of the interac- 
tion of the scalar field with the time-dependent potential vector. Here, the 
frequency is time-dependent because the scalar field lives in a time-dependent 
background, or, in some sense, because the scalar field interacts with the clas- 
sical gravitational background. Therefore, we already see at this stage that we 
can have particle creation due to the interaction with a classical gravitational 
field (instead of a classical electric field in the previous Section). Of course, 
the two cases are not exactly similar in the sense that the time dependence of 
uj 2 is different. Indeed, in the Schwinger case, uj 2 (k, t) typically contains terms 
proportional to t and t 2 , see Eq. ([8]), while, in the inflationary case, the term 
a" I a is typically proportional to 1/rj 2 . As a consequence, the solution to the 
mode equation and the particle creation rate will be different even if, again, 
the basic mechanism at play is exactly the same in both situations. 

The mode amplitude fi k is complex but one can also work with real vari- 
ables fi k and fx k , as was done previously in Eq. (I10|) . defined such that 

^ k = 71 + l ^ k) ' (48) 

In terms of these variables, the relation fi k = \L- k reads [i k = /j^ k and 
fi k = —fJ. I _ k - Then, the action (or Lagrangian) of the system takes the form 

s=± [d v [ d 3 fc{( M r) 2 + (ML') 2 - 2- + 

* J JR 3 + I a 

n' 2 

^-k 2 



(t4) +(^r . (49) 



One can check that it also leads to the correct equations of motion for the 
two real variables fi k and fi k . 
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We can now pass to the Hamiltonian formalism. The conjugate momentum 
to \i k is defined by the formula 

5C , a! . 
Pk = J^=^-^ k - (50) 

One can check that the definitions of the conjugate momenta in the real 
and Fourier spaces are consistent in the sense that they are related by the 
(expected) expression 

n (v,x) = -0j,f<lkp k e* x . (51) 

We see that the definition of the conjugate momentum p k as the derivative of 
the Lagrangian in Fourier space with respect to /i' k * and not to n' k is consistent 
with the expression of the momentum in real space. Otherwise the momentum 
n(rj, x) in real space would have been expressed in terms of p k instead of p k . 
Moreover, one can also check that 

5C „, a' „ , 8C .a' , _ 

Mg'-T^' ^ 7-77 = /4' - -A4 , (52) 



" fc " 8 nl' ™ a ^ k ' rk ~ 5 ft ™ a 
and, clearly, we have 

Pk = ^ (Pg + ip k ) , (53) 

as expected. 

We are now in a position where we can compute explicitly the Hamiltonian 
in the momentum space. The Hamiltonian density, TL, is defined in terms of 
the Hamiltonian H of the system through the relation 

H= f d 3 kH = [ d 3 k (p fe/ 4* +p* kl i' k - C) , (54) 

and we obtain 

W = PkP* k + k 2 p k n* k + — (Pfe/Xfe + Pk^k) ■ (55) 

Let us make some comments on this expression. If the background gravita- 
tional field is not time-dependent, that is to say if the scalar field lives in 
Minkowski space-time where a' = 0, then the above Hamiltonian reduces 
to a free Hamiltonian : there is simply no classical "pump field" . From the 
Schwinger effect point of view, this would be similar to a situation where there 
is no external classical electric field. In these two cases, no particle creation 
would occur. Moreover, one can also check that the Hamilton equations 

d V - dp k ~ CPk+ a^' dr, d Mfe ~ a Pk c Mfc ' (5&) 
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lead to the correct equation of motion given by Eq. (|46p . Finally, in terms of 
the real variables, the Hamiltonian density reads 



H 



to) + + ^ 2 W 



(57) 

We notice that 7i is simply the sum of two identical Hamiltonians for para- 
metric oscillator, one for ^ and the other for /i : fc . 

The expressions (|55|) and (|57|) should be compared to Eqs. ([9]) and (fTTjl . 
We see that, although similar, the formulae are not identical. However, as we 
are now going to show, this difference is only apparent. Indeed, let us now 
restart from the Lagrangian given by Eq. (|49|) . One can always add a total 
derivative without modifying the underlying theory. If one adds the following 
term 



2d V [a [ ^ k) 
then the Lagrangian takes the form 



a I I \2 



(58) 



d a HW + 0O -w 2 (M) 0|) + (Mfc) ,(59) 



where w 2 (fc,?y) = k 2 — a" /a. In this case, the conjugate momenta are simply 
given by p£ = /i^' and p l k — fi k '. As a consequence, the Hamiltonian now 
reads 



H = 



R3+ [ Z Z Z 



(60) 



This time, the Hamiltonian is exactly similar to the Schwinger Hamiltonian 
given by Eq. (fTTj) . This is another manifestation of the fact that, except for the 
exact time dependence of the effective frequency, the physical phenomenon, 
namely particle creation under the influence of an external classical field, is 
the same in both cases. 

Let us now investigate in more detail the relation between the Hamiltonian 
given by Eq. (|5T|) and the one of Eq. (f6"0"]) . We have just seen that the two 
corresponding theories differ by a total derivative and, hence, are physically 
equivalent. Another way to discuss the same property is through a canonical 
transformation. For this purpose, let us consider the following Hamiltonian 



Hi (pi,?i) = \^p\ 



1 



-p\q\ 



- 1 2 2 



(61) 



where 0(77) is an arbitrary function of the time. Clearly, H\ plays the role of 
the Hamiltonian in Eq. l[57j) and a{rj) is the scale factor. Then, let us consider 
a canonical transformation of type II [29) such that (<?i,pi) — > (92,^2)) the 
generating function of which is given by (a similar transformation has also 
been studied in Refs. [501 [5T] ) 
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G 2 (qi,P2,v) = 9iP2- ll- ( 62 ) 
I a 

From this function, it is easy to establish the relation between the "old" 
variables and the "new" ones. One obtains 

dG 2 a' dG 2 

Pi = -5 — =P2 qi , Q2 = -5 — = 9i ■ (63 

oqi a dp 2 

In particular, the first relation reproduces Eq. (|50|) with the correct sign. 
Finally, the "new" Hamiltonian is given by 

H 2 ( P2 ,q 2 ) = Hl + ^ = 1 -pt+ 1 -(k 2 - q\ . (64) 

Clearly this Hamiltonian is similar to the Hamiltonian of Eq. (|6"0|). 

Therefore, the two versions of the theory, the one given by the Hamil- 
tonian ([57]) . which is what we naturally obtain in the case of cosmological 
perturbations (see Sec. [4]), and the one which leads to the Hamiltonian ([60]) 
"a la Schwinger" are simply connected by a canonical transformation and, 
thus, are physically identical. In the following, we will see that this is also the 
case at the quantum level. 



3.2 Quantization and the squeezed states formalism 

So far, the discussion has been purely classical. We now study the quantization 
of the system starting with the Heisenberg picture. The quantization in the 
functional picture that we used for the Schwinger effect will be investigated 
in the next sub-Section. At the quantum level, [Lu and pk become operators 
satisfying the commutation relation 



[£fc,4] =^ (3) (fe-p) 



(65) 



Clearly, factor ordering is now important. The quantum Hamiltonian is ob- 
tained from the classical one by properly symmetrizing the expression (|55 j) . 
This leads to 



H 



[ d 3 fc p k p{ + k 2 fi k (l{ + %-C 



Pfe/4 



PfeAfc + Pkp{ 



(66) 



In addition, this guarantees the hermiticity of the Hamiltonian. The next 
step consists in introducing the normal variable dk [28] (which becomes the 
annihilation operator, cl becoming the creation operator) defined by 



Ck(v) 



Pk ■ 



(67) 



Equivalently, one can also express /tfc and Pk in terms of the normal variable 
and its hermitic conjugate. This gives the following two relations 
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J_(, 
V2k V 



Cfe + ci k 



Pk 



(68) 



Then, from the commutation relation (165|) . or equivalently from the relation 
in real space [S(rj, x), II(r], y)\ — iS^(x — y), it follows that [ck(v) i c p{v)] = 
5^ (k—p). In terms of the normal variables, the scalar field and its conjugate 
momentum can be expressed as 



(p(rj, k) 
n{r),x) 



a{ri) (2tt) 3 / 2 
(2tt) 3 / 2 



dfe 



dfc i 



2k: 
k 



c k (v)e lk x - c{( V )e-* k x 



(69) 
(70) 



Obviously, when a' = 0, we recover the flat space-time limit and so we expect 
the time dependence of the normal variables to be just Ck(rj) oc e lkv . 

We can now calculate the Hamiltonian operator in terms of the creation 
and annihilation operators. Using Eq. (|66p one obtains 



H = 



1 



d A k 



k (cfcCfc 



)-<-(■ 



CkC-k - C-k^Ck^) 



(71) 



where it is important to notice that the integral is now calculated in R 3 and 
not in R 3+ . Let us analyze this Hamiltonian. The first term is the standard 
one and represents a collection of harmonic oscillators. The most interesting 
part is the second term. This term is responsible for the quantum creation of 
particles in curved space-time. It can be viewed as an interacting term between 
the scalar field and the classical background. The coupling function ia! /a is 
proportional to the derivative of the scale factor and, therefore, vanishes in 
flat space-time. From the structure of the interacting term, i.e. in particular 
the product of two creation operators for the mode k and — fc, we can also see 
that we have creation of pairs of quanta with opposite momenta during the 
cosmological expansion (thus momentum is conserved as it should), exactly 
as we had particle creation due to the interaction of the scalar field with a 
classical electric field in the previous section. 

We can now calculate the time evolution of the quantum operators (here, 
we are working in the Heisenberg picture). Everything is known if we can 
determine the temporal behavior of the creation and annihilation operators; 
this behavior is given by the Heisenberg equations which read 



defe 



c k ,H 



drj 



Cfe', 



(72) 



Inserting the expression for the Hamiltonian derived above, we arrive at the 
equations 



.dcfc a 

l— — = KCfc + l — C- k 

dr] a 



.dcfc t t .a' 

l— = —KCk + I — C_ 

drj a 



k ■ 



(73) 
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This system of equations can be solved by means of a Bogoliubov transfor- 
mation and the solution can be written as 

Cfe(f?) = u k {rj)c k {r)im) + «fe(?7)c-fc t (?7ini) , (74) 
CkHv) = ut(v)ckHmm) + v%(r])c- k (r]- m i) , (75) 

where ?7i n i is a given initial time and where the functions Uk{rf) and v k (rj) 
satisfy the equations 

. dukjrj_ _ kUk ^ + i—v*^) , ^hiv) = k Vk (ri) + i—u* k {rj) . (76) 
aq a drj a 

In addition, these two functions must satisfy |itfc| 2 — \v k \ 2 = 1 such that 
the commutation relation between the creation and annihilation operators is 
preserved in time. A very important fact is that the initial values of u k and 
v k are fixed and, from the Bogoliubov transformation, read 

UkiVini) = 1 , Vfefeni) = . (77) 

Therefore, we remark that, in some sense, the initial conditions are fixed by 
the procedure of quantization. In fact, Eqs. (|77| imply that the initial state 
has been chosen to be the vacuum |0) at time r\ — r]^. A priori, it is not 
obvious that this choice is well-motivated but it turns out to be the case in 
an inflationary universe. This property constitutes one of the most important 
aspect of the inflationary scenario. Here, we do not discuss further this issue 
but we will come back to the problem of fixing the initial conditions at the 
beginning of inflation in the following. 

At this point, the next move is to establish the link between the formalism 
exposed above and the classical picture. For this purpose, it is interesting to 
establish the equation of motion obeyed by the function Uk + Straightfor- 
ward manipulations from Eqs. (|76j) lead to 

(«* + »£)" + (fc 2 -^) («* + <) = 0. (78) 

We see that the function u k + v k obeys the same equation as the variable /ik ■ 
This is to be expected since, using the Bogoliubov transformation, the scalar 
field operator can be re-written as 



A, - 1 1 r dk 



a(r,) (2tt) 3 / 2 J ^2k 



(u k +v* k ) (r/)ck(vim)e lk x 

(79) 



Therefore, if we are given a scale factor a(j]), we can now calculate the com- 
plete time evolution of the quantum scalar field by means of the formalism 
presented above. 
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In fact, the Bogoliubov transformation (|74p and ([75]) can be expressed in 
a different manner which is useful in order to introduce the squeezed states 
formalism. For this purpose, let us come back to the functions Uk and Vf.. We 
have seen that, in order for the commutator of the creation and annihilation 
operators to be preserved in time, these two functions must satisfy \uk\ 2 — 



1. This means that we can always write 



Uk — e l k cosh rk 



Vk 



= c -*( s x- 2 <t»<) S inh 



rk 



(80) 



where the quantities r^, 9k and <f>k are functions of time. They are called the 
squeezing parameter, rotation angle and squeezing angle respectively. These 
functions obey the equations 



dr fe a' 
— — = — cos 24> k 
arj a 



— — - = — k sin 2(f>k coth 2rk , 

arj a 



— — = — k sm Z(pk tann ■ 

arj a 



(81) 
(82) 



These expressions can be used for an explicit calculation of r^, 9k and <pk 
when a specific scale factor a(rj) is given. Now, the crucial property is that 
the Bogoliubov transformation (|Ti|) , ([75]) which solves the perturbed Einstein 
equations can be cast into the following form (32 [33l [3H [35] 



c k ( v ) = R(e)S(r, l p)c k (r h )SHr,i P )R^e), 

where the operators R{9) and S(r, ip) are given by 

R(9) =cxp|-i6» fc 4(?7i)c fe (?7i) +c 1 ^ fe (?7 i )c_ fc (7 ?i ) 
S(r,<p) = exp(r fc [e- 2 ^4(7 7i )c fc (ry i ) - e 21 ^ cl k ( Vi )c_ k (Vi) 



(83) 
(84) 



(85) 
(86) 



Eqs. ([83]) and ([84]) allows us to interpret the Bogoluibov transformation in a 
new manner: indeed we can also see the time evolution of the creation and 
annihilation operators as rotations in the Hilbert space. 

The previous considerations are valid in the Heisenberg picture. What 
happens in the Schrodinger picture where the operators no longer evolve but 
the states become time-dependent? For the sake of simplicity, let us ignore 9k 
and cj)k by setting 9k = 4>k — 0. As mentioned above, let us also postulate that 
the system is originally placed in the vacuum state |0). Then, the previous 
results imply that, after the cosmological evolution, the mode characterized 
by the wave-vector k will evolve into the following state [221 1331 1311 135] 



exp< r k 



(87) 
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which is, by definition, a two-mode vacuum squeezed state. This state is a 
very peculiar state and is of particular relevance in other branches of physics 
as well, most notably in quantum optics 36J. 

We now discuss the properties of such a quantum state. For this purpose, it 
is interesting to recall that a state containing a fixed number of particles, \n), 
can be obtained by successive action of the creation operator on the vacuum. 
Explicitly, one has 

(J \n 

|n> = ^r|0). 



Let us also introduce the coherent (Glauber) quantum state [55]. It is defined 
by the following expression 



oo 



£^=N, (89) 



where a is a complex number. The coherent state is especially important 
in quantum optics since they represent, in a sense to be specified, the most 
classical state. We will come back to this question in the last Section. Finally 
we also define two new operators B and T by 



Bu — — 



{c k f - (c fe t) 2 



T k = e M " , (90) 



where r is a real number (in fact, our squeezing parameter). The operators Bk 
and Tfc possess various interesting properties, in particular Bk is anti-unitary, 
B' k = —Bk, and, as a consequence, Tk is unitary, T^T^ = 1. The general 
definition of a squeezed state |s) is given by 



= Tl\a)=e- B >°\a). (91) 

Let us notice that this is the expression for a one-mode squeezed state while, in 
Eq. (|87p , we have to deal with a two-mode squeezed state (hence the presence 
of operators c& and C-k that arises from the fact that we have pair creation, 
while in the above definition we only have operators c 2 . and cj, 2 ). The proper- 
ties of one and two-mode squeezed states are similar and, here, for simplicity, 
we focus on the one-mode state only. Moreover, in our case, \a) = |0) which 
means that, in the cosmological case, we have a two-mode vacuum state. 

Why is this state called a squeezed state? To answer this question, we 
introduce two new operators that are linear combinations of the creation and 
annihilation operators, namely 

(cfc)p = ^ (c k + 4) (Cfe)<2 = Yi ( Cfe ~ c fc) ' ( 92 ) 



These new operators are annihilation and creation operators of standing waves 
since, in a Fourier expansion of the field, they would stand in front of cos krj 
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and sin k-q rather than % %kr> and e v in the case of the standard creation and 
annihilation operators. Then, it is straightforward to demonstrate that 

(s\(c k ) P \s) = ^-e 1 - , (s\(c k ) Q \s) = . (93) 

Let us now calculate the mean value of the squares of these operators. We 
have 

2r 

(s\(c k ) 2 P \s) = ^-(a 2 + a* 2 + 2aa* + 1) , (94) 

and a similar expression for (s|(cfe)g|s) (but with e~ 2r instead of e 2r ). We are 
now in a position where the dispersion in the squeezed state of the operators 
(cfe)p and (cfc)g can be calculated. One finds 

A{c k ) P = <J (s\(c k ) 2 P \s) - ( S \(c k ) P \ S )i = e - , A{c k ) Q = e — : (95) 

and, therefore, from these equations one deduces that 

A(c k ) P A{c k ) Q = i . (96) 

We see that the lower bound of the Heisenberg uncertainty relations is reached 
but, contrary to a coherent state, the dispersion is not equal for the two opera- 
tors. On the contrary, the dispersion can be very small on one component and 
very large on the other hence the name "squeezed state" . In the cosmological 
situation, this is actually the case. Indeed, Refs. 32, 33, 34] have shown that, 
for modes whose wavelengths are of the order of the Hubble length today, that 
is to say the modes that contribute the most to the "large angle" CMB mul- 
tipoles Gn (corresponding to a frequency of uj ~ 10~ 17 Hz), one has r ~ 120. 
From Eqs. (|95p . we see that this corresponds to a very strong squeezing, in 
fact much larger than what can be achieved in the laboratory [33]. 

It is also clear that a strongly squeezed state is not a classical state in the 
sense that it is very far from the coherent state for which A(c k )p — A(c k )Q. 
On the other hand, since the mean value of N k = c k c k is given by 

(s\N k \s) =sinh 2 r, (97) 

a strongly vacuum squeezed state contains a very large number of particles 
and this criterion is often taken as a criterion of classicality. Therefore, we 
see that the meaning of classicality for a strongly squeezed state is a subtle 
issue [UJ [T31 [22] [33] [331 133 13H] since different criterions seem to give different 
answers. We will come back to this point in the last Section of this review 
article. 

3.3 Quantization in the functional approach 

Let us now discuss the quantization in the functional approach where each 
Fourier mode is described by a wave-function. For this purpose, we use the 
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description in terms of real variables. This will allow us to emphasize again 
the complete analogy that exists between the Schwinger effect and the the- 
ory of inflationary cosmological perturbations of quantum-mechanical origin. 
We restart from Eq. (|57J) and, since we deal with quantum operators, we sym- 
metrize the corresponding expressions. In this case, the quantum Hamiltonian 
reads 



H= I d J fc 

3 + 



a' k 2 

Y a + pl kFk) + y ^ 2 



H a +H 1 . (98) 



We also have the following commutation relations that are compatible with 
Eq. (JBSJ) 

[Afe ,P$]= ^ (3) (fc - P) > [A 1 * . = w (3) (fc - P) ■ (99) 

In the Schrodinger picture, similarly to Eqs. (|13p . the above-mentioned oper- 
ators admit the following representation 

A^ = /^, ^ = -i^g. (ioo) 

Therefore, one deduces that the Hamiltonian (here, the Hamiltonian for the 
real part of /Ltfe, hence for a fixed Fourier mode) can be written as 

1 d 2 W i a' a' d\P k 2 o 

Again, if a' = 0, we recover the Hamiltonian of an harmonic oscillator (instead 
of the Hamiltonian of a parametric oscillator when a' ^ 0) . 

Let us now study the ground state of the theory. As done in Eq. ([TBI) , we 
have the following Gaussian state, 

^(Wt)=*(l)^' 4(,)W)! , (102) 

where Nh and are two functions to be determined. They are found by 
means of the Schrodinger equation id^u = HJ* &u that leads to 

N 1 i n' n' h 2 

# = ^ fc ^-> tf fc = -2i^-2-fl fc + i* (103) 
iVfc la a 2 

The analogy with Eqs. (fT7|) is obvious. We notice, however, that the structure 
of the equations is not exactly similar. This is due to the presence of the terms 
proportional to a' /a in the Hamiltonian (|101[) that have no equivalent in the 
Hamiltonian (fT5]) . Below, we briefly come back to this point. These equations 
can be integrated and the solutions read 
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%-(»*y\ a (104) 

\ t J 2 (/fe/aj 

where obeys the equation f k + (fc 2 — a"/a) = 0. Therefore, the inte- 
gration of the equation controlling the time evolution of the mode function 
leads to a complete determination of the quantum state of the system in full 
agreement with what was discussed before in the case of the Schwinger effect. 
Again, the fact that the solution for flk is given in terms of the function fk /a 
and not only in terms of fk, as one could have guessed from Eq. ([15]). is due 
to the presence of the terms proportional to a' /a in Eq. (|101|) . 

Let us briefly come back to the derivation of the above solution. Eq. (|103|) 
is a Ricatti equation and, therefore, can be solved in the usual way, namely 
by transforming this non-linear first order differential equation into a linear 
second order differential equation. In order to find Nk, one requires that the 
wave function is normalized, that is to say 

/Wd/*£ = i> ( 105 ) 

which leads to the previous expression of Nk- Moreover, there is also the 
following consistency check. The real part of the second of Eqs. (|103[) reads 

(SR/4)' = msik x 3X? fe - 2— sRf2 k , (106) 

a 

and the imaginary part of the first of Eqs. (|103[) can be written as N' k /Nk = 
^sQk — a'/ (2a). It is straightforward to check that, inserting the above solution 
for Nk into the last equation, precisely leads to Eq. (|106[) . 

Let us now come back to the remark made before that the structure of 
Eqs. (|103p is not exactly similar to what we have in the Schwinger case due to 
the presence of the terms proportional to a' /a. The reason is clearly that we 
have used the Hamiltonian given by Eq. (|57[) which contains such terms. But, 
obviously, one can also use the Hamiltonian given by Eq. (|60p. Then, assuming 
again the Gaussian form (| 1 02|) for the wave- function, the Schrodinger equation 
reduces to 

N' i 
i^ = ^k, f2' k = -2in 2 k + -u 2 (k,r 1 ). (107) 

which are now exactly similar to Eqs. (TP?)) . As a consequence, the solutions 
are also the same and read 

/29W2 fc \ 1/4 if' . 

»> = {-^r) . ^ = -27*- (108) 

where fk obeys the mode function equation f' k ' + cu 2 fk = 0. 

In Sec. 13.11 we have established, at the classical level, the equivalence 
between the two formulations discussed above, that is to say the one based on 
the Hamiltonian (|57|) . which leads to a Gaussian wave- function with Nk and 
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Qk given by Eqs. (|104|) . and the one based on the Hamiltonian (|60p. which 
also leads to a Gaussian wave-function but with and i?*, now given by 
Eqs. (|108p . We now study this link at the quantum level and, for this purpose, 
we reconsider the simple model introduced after Eq. (|61[) . In Sec. 13.11 we 
showed that the two formulations are connected by a canonical transformation 
and the question is now to implement this canonical transformation at the 
quantum level [301 EH EH [40], [41] H2]- For this purpose, one must find a 
unitary operator IA such that the relations 



'12 



Uqxrt , f>2 = Upirt , (109) 



exactly reproduce the classical analogues ([D3]) . A natural candidate would be 
the following operator 



U = e lG2 = exp 



- — Qi + ~ [QiPi + P1Q1) 
2 a I 



(110) 



where G 2 is generating function introduced in Eq. (|62p . However, as already 
remarked in Ref. [30], this choice is too naive and does not work. In order to 
understand what is going on, let us introduce a generalized version of Eq. (f6l"j) , 
following Eq. (2.21) of Ref. [3D], which at the classical level reads 

Hi (pi, qi) = \fov\ + 02(v)PiQi + \pik\l , (HI) 

where for simplicity we consider that [3^ and (3\ are constant while (3 2 is a time- 
dependent function (in Ref. [30] , all the 0i's are time-dependent functions). 
Clearly, our case corresponds to 0y = 0s = 1 and ,$2 = a' /a. Then, as before, 
we consider a canonical transformation of type II such that (qi,pi) — * {q 2 ,p 2 ) 
with the following generating function 

G 2 (quP2, rj) = ^ l! \iP2 - ^q\ ■ (112) 

2/J3 

Setting /3i = /?3 = 1 and /?2 = a'/ a i n the above expression reproduces 
Eq. (|62p as expected. Performing standard calculations, one finds that the 
relation between the "old" variables and the "new" ones reads 

dG 2 n-xi-i dG 2 n -l/2 

Pi = ^ — = P3 P2- -5-91 , <?2 = p — = 03 qi, (113) 
oqi 03 op 2 

and that the "new" Hamiltonian can now be expressed as 

H 2 (p 2 ,q 2 ) =H l + ^ = \pl + \ (ftftfc 2 - 0' 2 - 01) q\ . (114) 

Notice, in particular, that the coefficient P3 is no longer present in the term 
pl/2. Then, in agreement with Eqs. (2.22) and (2.23) of Ref. [3D], let us 
consider the following operator 
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U(qi,Pi,v) = exp 



' 2^ 2 9i ) ex P 



•jMi) 



(115) 



Inserting this operator in Eqs. (|109p and using the Baker-Campbell-HausdorfT 
formula, e A Be~ A = B + [A, B] -\ , leads to the transformation 



Pi 



- /3 3 ~ 1/2 



/?2 * * fl- 1 / 2 - 

P3 



(116) 



namely exactly Eqs. (|113p . but now at the quantum level. Therefore, we 
conclude that U in Eq. (|115|) is the operator generating the correct quan- 
tum canonical transformation. In addition, as one can check with the help of 
Eq. (1112|) , this operator is different from e lG2 , in particular due to the pres- 
ence of the factor In 0$. Let us also notice that a similar operator has been 
considered recently in Refs. [3TJGHI], which carries out an investigation very 
relevant for what is discussed here, and that a factor akin to ln/?3 was also 
present in the operator U of that paper [see Eq. (2.46) where this factor is 
written as "ln\/l2/a"]. Moreover, and this is the main reason why we have 
considered a generalized version of Eq. (|6ip , we notice that our case is in fact 
very special since it corresponds to fa = 1 or In 03 = (or "e = 0" in the 
language of Ref. [39]). This means that, in the operator (|115|) . the second ex- 
ponential totally "disappears" while, of course, the term proportional to qiPi 
remains present in the classical generating function. Therefore, in our case, 
the quantum generating operator is just given by 



U(qi,Pi,r]) = exp ( ~-—qf 
Z a 



(117) 



Clearly, one can repeat the above calculations using this operator and show 
that this leads to Eqs. (|6"3")l but at the quantum level. 

Let us now turn to the transformation of the wave-function itself. It is 
given by 

f ((ft) = N 2 e~ n ^ =U\q ll p 1 ,r 1 )V(q 1 )=U\q 1 ,p 1 ,7 1 )N 1 e- n ^ , (118) 
from which, using Eq. (|117|) . one deduces that 



N 2 = m, Qi = Qi- i- 



(119) 



The relation N 2 = N-y comes from the fact that the quantity is given by 
the real part of the function (2\ 2 and that fl 2 and fl\ differ by a complex factor 
only. The above relation exactly reproduces what was observed in Eqs. (|104p 
and Eqs. (|108|) . From these formulae, we see that Q\ — —i/2(f /a)' /(/ /a), 
see Eqs. (|104j) . while Q 2 = —if/(2f), see Eqs. (|108[) . and they indeed satisfy 
Eq. (|119p . Of course, the wave- functions after the quantum canonical trans- 
formation is normalized because = UW^l^i = the operator U 
being unitary. 



26 Jerome Martin 

3.4 The Power Spectrum 

Let us now calculate the two-point correlation function in the quantum state 
where the scalar field is put by the cosmological evolution. As will be discussed 
in the following, this quantity is relevant in astrophysics because, in the case of 
cosmological perturbations, it is directly observable; in particular it is directly 
linked to CMB fluctuations. Its definition reads 

/n 

(120) 

In this formula, the brackets mean the quantum average according to the stan- 
dard definition, i.e. (A) = J dx <?* A{x) Then, using the Fourier expansion 
of the scalar field and permuting the integrals, one obtains 



S(r], x)&(rj, x + r) \ = 



1 1 

t? (2^)3 



dpdqe' 1 ' ~e 



ip x iq-(x+r) 



£r /23W2 fc 



n 



V n 



(121) 

The above expression vanishes unless p = —q. Indeed, if \p\ ^ \q\ then the 
quantity [i p [i q is "linear" in [ip 1 and fig' 1 and, consequently, the Gaussian 

integral is zero. If p = q, then /J, p /J, q = (^p) 2 - (m p ) 2 an d cacn term is indeed 
non- vanishing but the sum is zero because of the minus sign. Therefore, one 
obtains 



S( V ,x)S(r,,x + r)) = / dpe^JI 



nd^d ML e- 2 ^^^[K) 2+ K) 2 ] te)\i22) 



The overall factor of 2 originates from the fact that the integral over /i k is 
equal to the integral over fi k . The next step is to perform the path integral. In 
the above infinite product of integrals, all of them are of the form "J dxe~ ax " 
except the one over ^ which is of the form "f dxx 2 c~ ax2 " . Using standard 
results for Gaussian integrals, one gets 



<P(i], x)S(i], x + r) 



= 1J- [ 

a 2 (2tt) 3 J 



dpc 



tp-r 



n(^)n 



vmn~ k 



n-l , /— 



(123) 
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The infinite product "Hfc " means a product over all the wave- vectors 
but p. Clearly, one can complete this product by inserting an extra factor 
\pnj ^/25R!7 p coming from the last term in the integral. Then, the last two 
products exactly cancel the first one. Finally, one obtains the simple expres- 
sion 

- hw? I dpelpr vk- (124) 



Using the form of i? p in the ground state wave- function, see Eq. (|108[) . one 
obtains 

2^ p = -i ^-^ = -JL-, (125) 

where we have used the fact that, with the initial condition (we will return 
to this point in the following) p p (rf) — > (2p)~ 1 ' 2 e ipn when prj — > — oo, the 
Wronskian is equal to i. 

At this point, one can also make the following remark about the canonical 
transformation discussed in the previous subsection. It is clear that the power 
spectrum must be the same before and after the canonical transformation. 
Above, we used the form of fi p given by Eq. (| 1 08[) . But one could have used 
the form given by Eq. (fT04|) in the same manner and without affecting the 
final result. Technically, this can be seen in Eq. (|124p where it is clear that the 
power spectrum only depends on $tf2 p . Since we demonstrated before that the 
canonical transformation only modifies the imaginary part of i? p , the power 
spectrum remains indeed the same. 

Therefore, the final expression reads 

$(r,,x)$(r,,x + r)) = J>^3 / dpe 4p V pMp (126) 

(127) 



1 f dp sin pr ^ 2 p p 2 



4ir 2 J p pr 



-v 2 



This expression is the standard one, usually derived in the Heisenberg pic- 
ture |16[ll7j. Knowledge of the mode function (including the initial conditions) 
is sufficient to estimate the power spectrum. In the following, we consider the 
case of inflationary cosmological perturbations and investigate which quan- 
tity plays the role of pk in that framework. This will allow us to discuss the 
inflationary predictions. 



4 Inflationary Cosmological Perturbations of 
Quantum-Mechanical Origin 

4.1 General Formalism 

In this section, we finally consider our main subject, namely the theory of 
inflationary cosmological perturbations of quantum- mechanical origin |11[ I16[ 
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[T71 33] . Our goal is to go beyond the isotropic and homogeneous FLRW Uni- 
verse, the metric of which can be written as 



ds 2 = a 2 (r]) 



-drf 



sifdx i dx j 



(128) 



and to study how small quantum perturbations around the above-mentioned 
solution behave during inflation. As we will see, the basic physical phenomenon 
and, hence, the corresponding formalism are similar to what was discussed 
before. As already emphasized, we are mainly concerned with inflation, that 
is to say a phase of accelerated expansion that took place in the early universe. 
In general relativity, such a phase can be obtained if the matter content is 
dominated by a fluid whose pressure is negative. Since, at very high energies, 
quantum field theory is the natural candidate to describe matter, it is natural 
and simple to postulate that a scalar field (the "inflaton" ) was responsible for 
the evolution of the universe in this regime. Therefore, the total action of the 
system is given by 



S 



m 

pi 
I6n 



d 4 a 



-gR 



d 4 3 



g^d^ip + V[fp) 



(129) 



where tp is the inflaton field. Our discussion will be (almost) independent of 
the detailed shape of the potential V(ip) but, clearly, deriving from high energy 
physics (for instance string theory) what this shape could be (in particular 
explaining the required flatness of the potential) is a major issue |46[ 147] . 

Beyond homogeneity and isotropy, the most general form of the perturbed 
line element can be expressed as [llj : 



ds' 



a 2 (?7){-(l - 2cf>)drf + 2(d l B)dx i drj + [(1 - 2ip)8, 



(•'!) 



hij\dx' l dx 3 } 



(130) 



In the above expression, the functions <f>, B, ip and E represent the scalar 
sector whereas the tensor hij, satisfying hi 1 = = hij ,J , represents the gravi- 
tational waves. These functions must be small in comparison to one in order 
for the perturbative treatment to be valid. There are no vector perturbations 
because a single scalar field cannot seed rotational perturbations. At the linear 
level, the two types of perturbations decouple and, therefore, can be treated 
separately. 

In the case of scalar perturbations of the geometry evoked above, the 
four functions are in fact redundant (thanks to our freedom to choose the 
coordinate system) and, in fact, the scalar fluctuations of the geometry can 
be characterized by a single quantity, namely the gauge-invariant Bardeen 
potential ^ B [43] (not to be confused with the scalar field <P considered before) 
defined by 



<P B (r,,x) = ^+-[a (B-E')}' 



(131) 
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On the other hand, the fluctuations in the inflaton scalar field are character- 
ized by the following gauge-invariant quantity 8(p & > 

8(pW fa, x) = 6<p + ip' (B - E') . (132) 

We have therefore two gauge-invariant quantities but only one degree of free- 
dom since <P B and Stp^ are coupled through the perturbed Einstein equa- 
tions. As a consequence, in the scalar sector of the theory, everything can 
be reduced to the study of a single gauge-invariant variable (the so-called 
Mukhanov-Sasaki variable) defined by [3] 



v(r),x) 



n 



(133) 



Let us notice that we will also work with the rescaled variable fi g defined by 
H s (77, x) = —\[2kv. Finally, density perturbations are also often characterized 
by the so-called conserved quantity C (v> x ) [SI US] defined by fj, s = —2ay^(, 
where 7 = 1 - H'/H 2 . 

In the tensor sector (which is automatically gauge invariant), the quan- 
tity which plays the role of fj, s (r),x) is /i T (77, x), defined according to hij — 
(fi T /a)Qij, where Qij are the (transverse and traceless) eigentensors of the 
Laplace operator on the space-like sections |43) . 

As usual, it is more convenient to study the perturbations mode by mode 
and, for this purpose, we will follow the evolution of the perturbations in 
Fourier space. Therefore, the study of cosmological perturbations during in- 
flation reduces to investigating the behaviors of only two variables: fi s k (77) 
and n T k (r))- 

Let us now establish the equations of motion for our two basic quantities. 
Since we want the variation of the action (|129p to give the first order equations 
of motion for fj, s k (77) and /x T fe (77), we have to expand the action pertubatively 
up to second order in the metric perturbations and in the scalar field fluctu- 
ations. After a lengthy and tedious calculation, one obtains [11] 



^6S = ~ / d 4 * 



(v'Y -5 lJ divdjV + 



9 

m 



/ d 4 xa 2 (77) [(h^)'(h^y - dkiWj^ihh)] , (134) 

Notice that the constant m pl does not appear explicitly in the scalar part of the 
action because it has been absorbed via the background Einstein equations 
(however, see also Ref. [21] )■ It is also important to stress again that the 
previous expression is valid for any potential V(<p). 

Variation of the action leads to the following equation of motion for the 
two quantities /i s k (77) and /_t T fc (77) 



d.77 2 



+ w s, T fc > J 7)Ms,Tfc = °> ( 135 ) 
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with 



u, s 2 (M) = fc 2 



(«V7) 
a^7 



(136) 



We have thus reached our goal and demonstrated that cosmological perturba- 
tions obey exactly the same type of equation as a scalar field interacting with 
a classical electric field (Schwinger effect), namely the equation of a paramet- 
ric oscillator as can be checked by comparing Eq. (|135[) with Eq. ([7]). The 
only difference lies in the physical nature of the classical source. In the case of 
cosmological perturbations, the (background) gravitational field is the classi- 
cal source. The time dependence of the frequencies co s and uj t is also different 
(recall that, in the case of the Schwinger effect, uj 2 contains terms proportional 
to t and t 2 ). Here, the dependence is fixed as soon as the behavior of the scale 
factor a(r/) is known. It is also interesting to notice that, a priori, the time 
dependence of u s is not the same as the one of lu t . Indeed, lu t depends on 
a and its derivatives up to second order while u s depends on the scale factor 
and its derivatives up to the fourth order (since it contains a term 7", the 
quantity 7 containing itself a term a"). Finally, the quantization of the theory 
proceeds as before and, as a consequence of the interaction between the quan- 
tum cosmological perturbations and the classical background, this results in 
the phenomenon of particle creation, here graviton creation. Classically, this 
corresponds to the amplification ("growing mode") of the fluctuations. 

In the next section, we describe this phenomenon for an inflationary scale 
factor. 

4.2 The Inflationary Effective Frequencies 

So far, we have never specified a(rj) and, a priori, the mechanism of graviton 
creation is valid for any scale factor provided it is time-dependent. However, 
clearly, the detailed properties of the transition amplitude (tfql^) depend on 
the time behavior of the effective frequency uj 2 {k,rf) and, hence, on the form 
of a(rj). Obviously, in the case of the Schwinger effect, a frequency different 
from the one given by Eq. ([8]) would have led to a number of created pairs 
different from Eq. ([39]) . 

In order to evaluate oj 2 {k, rj) and u> 2 (k, rj) for a typical inflationary model, 
one can use the slow-roll approximation E01 [S3 HI] ■ Indeed, during infla- 
tion and by definition, the kinetic energy to potential energy ratio and the 
scalar field acceleration to the scalar field velocity ratio are small and this 
suggests to view these quantities as parameters in which a systematic expan- 
sion can be performed. Therefore, one introduces the two parameters e\ and 
€2 [22] according to 



From the above expressions, one sees that ei/3 measures the ratio of the 
kinetic energy to the total energy while £2 > (respectively €2 < 0) represents 



= 3 




(137) 
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a model where the kinetic energy itself increases (respectively decreases) with 
respect to the total energy. It is also interesting to notice that e 2 = 2ei marks 
the frontier between models where the kinetic energy increases (e 2 > 2ei) and 
the models where it decreases (t2 < 2ei). Provided the slow-roll conditions 
are valid, that is to say ex,2 1, one can also invert the previous expressions 
and express the slow-roll parameters only in terms of the inflaton potential. 
This leads to 



ei 



16tt 



C-2 



47T 



v 



(138) 



where, in the present context, a prime denotes a derivative with respect to 
the scalar field (p. Concrete calculations of slow-roll parameters for specific 
models can be found in Ref. [16] . 

Then, one can show that the two effective frequencies, to first order in the 
slow-roll parameters, can be expressed as [TTJ 05] 



, 2 (k, V )~k 2 - 



2 + 3e x - 3e 2 /2 



W 2(M)-fc 2 - 



2 + 3ei 
V 2 



(139) 



Several remarks are in order at this point. Firstly, and as already mentioned 
previously, the time dependence in the inflationary case is different from the 
Schwinger case: the effective frequency contains terms proportional to 1/r] 2 . 
Therefore, although the basic physical phenomenon is the same, one can ex- 
pect the detailed predictions to differ. Secondly, different inflationary models 
correspond to different inflaton potentials (or to different time variations of 
the scale factor) and, hence, to different values for the slow-roll parameters. 
One notices that the effective frequencies are sensitive to the details of the 
inflationary models since oj 2 {k,rf) and w^(fc,?7) depend on t\ and £2. 



4.3 The WKB Approximation 

We have established the form of the effective frequencies in the case of infla- 
tion. One must now solve the mode equations (|135[) . For this purpose, we now 
reiterate the analysis of Sec. 12.31 using the WKB approximation [53]. As was 
the case for the Schwinger effect, the mode function can be found exactly. It 
is given in terms of Bessel functions [instead of parabolic cylinder functions, 
see Eq. flU] 

/ifc('7) = Vkv [A h Ju {krj) + B k J-„ {kij)} , (140) 

where the orders are now functions of the slow-roll parameters, v s = —3/2 — 
ei — e 2 /2 and v T = — 3/2 — ei. Then, one must choose the initial conditions. As 
discussed in the case of the Schwinger effect, we use the WKB approximation 
to discuss this question. The first step is to calculate the quantity Q in order 
to identify the regime where an adiabatic vacuum is available. Straightforward 
calculations lead to 
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Fig. 2. Evolution of the quantity \Q/u 2 \ with the quantity kr] for a typical model of 
inflation according to Eq. (|141|l (we have neglected the corrections proportional to 
the slow-roll parameters). In the limit kr] — *• -co, which corresponds to a wavelength 
much smaller than the Hubble radius, |Q/u; 2 | vanishes and the notion of an adiabatic 
vacuum is available. 
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UJ 2 
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1 - 3k z r] 



2„,2 



(l-fcV/2) 



(141) 



This quantity is represented in Fig. [2] 

Let us now discuss this plot in more detail. The problem is characterized 
by two scales: the wavelength of the corresponding Fourier mode given by 



(142) 



where k is the co- moving wavenumber, and the Hubble radius which can be 
expressed as 



^ (v) 



(143) 



We notice that the quantity |Q/ct> 2 | vanishes in the limit kr] — > — oo. This 
limit corresponds to a regime where A <C £ n ■ In this case, the wavelength is so 
small in comparison with the scale £ a characterizing the curvature of space- 
time that the Fourier mode does not feel the FLRW Universe but behaves as 
if it were in flat (Minkowski) space-time. Clearly, in this regime, an adiabatic 
vacuum state is available since we recover the standard quantum field theory 
description. In the limit kr] — > 0, the quantity |Q/w 2 | goes to 1/8 = 0(1) as 
can be checked in Fig. O This regime corresponds to the case where A > l H , 
that is to say when the wavelength of the Fourier mode is outside the Hubble 
radius. In this case, the curvature of space-time is felt and, as a consequence, 



The Cosmological Schwinger Effect 



33 



the WKB approximation is violated and there is no unique vacuum state in 
this limit. 

We have just seen that when a mode is sub-Hubble, that is to say A -C £ n , 
the WKB approximation is valid. Let us notice that, without a phase of in- 
flation, all the Fourier modes of astrophysical interest today would have been 
outside the Hubble radius in the early Universe. It is only because, during 
inflation, the Hubble radius is constant that, initially, the Fourier modes are 
inside the Hubble radius. Therefore, although it was not designed for this 
purpose, a phase of inflation automatically implies that the WKB approxima- 
tion is valid in the early Universe and, as a consequence, ensures that we can 
choose a well-defined initial state. This is clearly an "extra bonus" of utmost 
importance. In the adiabatic regime, the solution for the mode function can 
be written as 

MfeO?) = UkUwkbAv) + /3fcMwkb,fc(?7) > ( 144 ) 

where 

= .„ \ e-^-"t fc ^. (145) 

VMM) 



As done for the Schwinger effect, see Eq. ([26]). we now choose the initial 
conditions such that au = 1, /3ft = 0, corresponding to only one WKB branch 
in Eq. (I144j) . This completely fixes the coefficients Ak and B k in Eq. ()140|1 . 
One obtains [compare with Eqs. ([35]) and (|36]) ] 



£ = -e*~ 5 fc = - — G - ?r -— . (146) 

Equipped with the above exact solution for the mode function, the inflationary 
predictions can be determined. 

Before turning to this calculation, let us quickly come back to the fact 
that the WKB approximation breaks down on super-Hubble scales. In fact, 
this problem bears a close resemblance with a situation discussed by atomic 
physicists at the time quantum mechanics was born. The subject debated 
was the application of the WKB approximation to the motion in a central 
field and, more specifically, how the Balmer formula for the energy levels of 
hydrogenic atoms, can be recovered within the WKB approximation. The 
effective frequency for hydrogenic atoms is given by (obviously, in the atomic 
physics context, the wave equation is not a differential equation with respect 
to time but to the radial coordinate r) 

" 2 ( E > = w [E + — j - A_J , (147) 

where Ze is the (attractive) central charge and £ the quantum number of 
angular momentum. The symbol E denotes the energy of the particle and 
is negative in the case of a bound state. Apart from the term Ze 2 jr and 
up to the identification r <-> 77, the effective frequency has exactly the same 
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form as oj s T (k, rj) during inflation, see Eqs. (|136[) . Therefore, calculating the 
evolution of cosmological perturbations on super-Hubble scales, |fcry| — > 0, is 
similar to determining the behavior of the hydrogen atom wave function in 
the vicinity of the nucleus, namely r — > 0. The calculation of the energy levels 
by means of the WKB approximation was first addressed by Kramers [54] 
and by Young and Uhlenbeck [55]. They noticed that the Balmer formula 
was not properly recovered but did not realize that this was due to a misuse 
of the WKB approximation. In 1937 the problem was considered again by 
Langer [53]. In a remarkable article, he showed that the WKB approximation 
breaks down at small r, for an effective frequency given by Eq. (fl47|) and, m 
addition, he suggested a method to circumvent this difficulty. Recently, this 
method has been applied to the calculation of the cosmological perturbations 
in Refs. [53] [57]. This gives rise to a new method of approximation, different 
from the more traditional slow-roll approximation. 

4.4 The Inflationary Power Spectra 

In this sub-Section we turn to the calculation of the inflationary observables. 
The first step is to quantize the system. Obviously, this proceeds exactly as 
for the Schwinger effect or for a scalar field in curved space-time, the two cases 
that we have discussed before. We do not repeat the formalism here. As before, 
in the functional Schrodinger picture, the wave-function of the perturbations 
is given by 

n n 

* r =n^(^^D=n^ R (^)^(^) > (148) 

k k 

with 

n fa.Mg) = N k (rj)e~ n ^(^) 2 , S$ ( V ,^ k ) = N k (77) e - n >W(»lY , 

(149) 

where the functions Nk(rj) and fikiv) are functions that can be determined 
using the Schrodinger equation. This leads to expressions similar to Eqs. (fl"7|) 
and (I107|) where now uj 2 (k, t) should be replaced by uj 2 t according to whether 
one considers the scalar perturbations or the gravitational waves. In particu- 
lar, the function f^kiv) is still given by — i// fc //ifc, see Eq. (|108[) . where, in the 
present context, ^ is given by the Bessel function of Eq. (|140p . 

At this stage, one could compute the amplitude (0~|0 + ) as one did in 
the case of the Schwinger effect. However, in the context of inflation, this 
is not the observable one is interested in. Indeed, we want to evaluate the 
amplitude of the fluctuations at the end of inflation and on super-Hubble 
scales. In this regime, as discussed before, there is no adiabatic state. So, in 
the context of inflation, there exists a "in" vacuum state |0 _ ) when krj — ► — oo 
but there is no "out" region and, consequently, no |0 + ) state. Of course, if 
one follows the evolution of the mode after inflation, then the unicity of the 
choice of the vacuum state is restored when the mode re-enters the Hubble 
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radius either during the radiation or matter dominated eras. But our goal is 
to compute the spectrum at the end of inflation. In other words, and contrary 
to the Schwinger effect, the quantity (0~|0 + ) is not really relevant for the 
inflationary cosmological perturbations. 

In fact, our goal is to calculate the anisotropies in the CMB (and/or to 
understand the distribution of galaxies) . The key point is that the presence of 
cosmological perturbations causes anisotropies in the CMB: this is the Sachs- 
Wolfe effect [58, 59] . More precisely, on large scales, one has 

^(e)ocC = -^, (150) 

where e represents a direction in the sky. The exact link is more complicated 
and has been discussed in details for instance in Refs. [171 [59]. I n fact, it is 
convenient to expand this operator on the celestial sphere, i.e. on the basis of 
spherical harmonics 

f~ +oo m— I 

^( e )=E E *tmY tm (9,tp). (151) 

e=2m=-l 

This allows us to calculate the vacuum two-point correlation function of tem- 
perature fluctuations. One gets 

(f (ex)f (e 2 )\ =£^C^(c 0S7 ) , (152) 

where Pf is a Legendre polynomial and 7 is the angle between the two vectors 
ex and e 2 . In the above expression, the brackets mean the standard quantum 
average. In practice, the observable two-point correlation function is rather 
defined by a spatial average over the celestial sphere. These two averages are 
of course not identical and the difference between them is at the origin of the 
concept of "cosmic variance" ; see Ref. [60] for a detailed explanation. The Ci 
's are the multipolc moments and have been measured with great accuracy 
by the WMAP experiment [5] . Clearly, as can be seen in Eq. (I150p , the above 
correlation function is related to the two-point correlation function of the 
cosmological fluctuations. Therefore, the relevant quantities to characterize 
the inflationary perturbations of quantum-mechanical origin are 

C( V , x)C(r,, x + r))^^ 00 f ^k 3 P c , (153) 

for scalar perturbations and, for tensor perturbations 

r+0 ° dk sinfcr. 



(h ij {r,,x)&(r,,x + r)) = J ^ h ^-k 3 P h , (154) 
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One can then repeat the calculation done in Sec. 13.41 in order to evaluate the 
above quantities. Indeed, the calculation proceeds exactly in the same way 
since the wave-functional is still a Gaussian. This gives 



k 3 P ( (k) 



8^2 



Ms 



k 3 P h (k) = 



M-i 



(155) 



These expressions should be compared with Eq. (|127p . 

The two power spectra can be easily computed using the exact solution for 
the mode function, see Eq. (|140p . At first order in the slow-roll parameters, 



one arrives at [49j [50l [51 
H 2 
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(156) 



(157) 



where C is a numerical constant, C ~ —0.73, and fc* an arbitrary scale called 
the "pivot scale". We see that the amplitude of the scalar power spectrum 
is given by a scale-invariant piece (that is to say which does not depend on 
k), H 2 / (ireim 2 ^ , plus logarithmic corrections the amplitude of which is con- 
trolled by the slow-roll parameters, namely by the micro-physics of inflation. 
The above remarks are also valid for tensor perturbations. The ratio of tensor 
over scalar is just given by k 3 Ph/k 3 P^ = 16ei. This means that the gravita- 
tional waves are always sub-dominant and that, when we measure the CMB 
anisotropies, we essentially see the scalar modes. This is rather unfortunate 
because this implies that one cannot measure the energy scale of inflation 
since the amplitude of the scalar power spectrum also depends on the slow- 
roll parameter t\. Only an independent measure of the gravitational waves 
contribution could allow us to break this degeneracy. On the other hand, the 
spectral indices are given by 



1 



In k 3 P, 



d\nk 



-2ei - e 2 



In k 3 Ph 



k=k. 



din A; 



-2ei. (158) 



k=k. 



As expected, the power spectra are always close to scale invariance (n s = 1 
and n T = 0) and the deviation from it is controlled by the magnitude of the 
two slow-roll parameters. 

To conclude this section, let us signal that the slow-roll parameters t\ 
and £2 are already constrained by the astrophysical data, see Fig. [3] for the 
constraints coming from the WMAP data. A complete analysis can be found 
in Refs. EHE51. 



5 The Classical Limit of Quantum Perturbations 



As discussed at length previously, the inflationary cosmological perturbations 
are of quantum-mechanical origin. However, from the observational point of 
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Fig. 3. 68% and 95% confidence intervals of the two-dimensional marginalized pos- 
teriors in the slow-roll parameters plane, obtained at leading order in the slow-roll 
expansion [5[Tj . The shading is the mean likelihood and the left plot is derived under 
an uniform prior on ei while the right panel corresponds to an uniform prior on 
log ei . 

view, it seems that we deal with a physical phenomenon where quantum me- 
chanics does not play a crucial role (even does not a play a role at all). 
Therefore, from the conceptual point of view, it is important to understand 
how the system can become classical [HI [T31 (see also Ref. [M] ) . We now 
turn to this question. 




5.1 Coherent States 

It seems natural to postulate that a quantum system behaves classically when 
it is placed in a state such that it follows (exactly or, at least, approximatively) 
the classical trajectory. For the sake of illustration, let us consider a simple 
one-dimensional system characterized by the Hamiltonian 

H(p,q) = ±p 2 + V(q), (159) 

where, for the moment, the potential V(q) is arbitrary. Solving the classical 
Hamilton's equations (given some initial conditions) 

provides the classical solution p c \ and g c j. At the technical level, the above- 
mentioned criterion of classicality amounts to choosing a state \W) such that 

Pd(t) = W(t)|#> , ?d(t) = • (161) 

This is clearly a non-trivial requirement as can be understood from the Ehren- 
fest theorem. Indeed, this theorem shows that, for any state one has 
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-{*\p(t)\&) = - (# 



OV(q) 



dq 



(¥\p\iP) , (162) 



These equations resemble the Hamilton's equations (|160j) but are of course 
not identical. This implies that, placed in an arbitrary state, the quantum 
system does not behave classically (i.e. the means of the position and of the 
momentum do not obey the classical equations). It would be the case only for 
a state |<P) such that 



dV(q) 



dq 



(163) 



which is not true in general since (fl/\q n \ty) ^ (\P\q\\P) n , but obviously satisfied 
if the potential assumes the particular shape V(q) oc q 2 , i.e. for the harmonic 
oscillator. In this case, the means of the position and of the momentum do 
follow the classical trajectory whatever the state I?/') is. This means that 
Eqs. (|16ip are in fact not sufficient to define classicality and that one needs 
to provide extra conditions. It seems natural to require that the wave packet 
is equally localized in coordinate and momentun Aq = y (q 2 ) — (q) 2 = Ap 
to the minimum allowed by the Heisenberg bound, namely AqAp = 1/2. 
This is another way to define a coherent state, see Eq. ([55)) which, therefore, 
represents the "most classical" state of a quantum harmonic oscillator. 

We now demonstrate that the state ([89]) indeed satisfies the above- 
mentioned properties. If the potential is given by V(q) = k 2 q 2 /2, then the 



Hamilton's equations can be expressed as p = 
mal variable" variable, see also Eq. (f6"T|) . 



-k q and q = p and the "nor- 



2 \ q+ k P 



(164) 



obeys the equation a — —ika which allows us to write the classical trajectory 
in phase space as 



Pd(t) = 

?ol(*) = 



with 



a 



ikt 



(a e~ 



■/ k I 



a*e lkt ) 



ikt\ 



q(t = 0) + -p(t = 0) 



(165) 
(166) 

(167) 



Let us consider that, at time t = 0, the system is placed in the state |«o)- At 
time t > 0, the integration of the Schrodinger equation leads to 



-ikt/2 



a {t) = a Q e- lkt ) 



(168) 
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This result should be understood as follows. In the expression ([89)1 which 
defines a coherent state, the factor a should be replaced with a e~ tkt to get the 
formula expressing the above state As already mentioned, at the quantum 
level, the normal variable becomes the annihilation operator [obtained from 
Eq. (|164p by simply replacing q and p with their quantum counter-parts] . This 
implies 

1 « f ) p = -i\ ~(a-al) . (169) 



2k y ' V 2 

Then, the crucial step is that any coherent state \a) is the eigenvector of a 
with the eigenvalue a, a\a) = a\a). Using this property, it is easy to show 
that, for the state |<?) defined by Eq. (|168|) , one has 

(W) = "<y| («oe-* fct - a* e ikt ) , (170) 
(*\q\*) = ^(<xoe- ikt + a* Q e ikt ) . (171) 



In the same way, straightforward manipulations lead to 

<f|p 2 |^) = ^{l + 43 2 [«W]} , (f|g 2 |fH^{l + 45R 2 [«(*)]} (172) 
from which one deduces 

We have thus reached our goal, i.e. we have shown that the state (j 168[) follows 
the classical trajectory and that the quantum dispersion around this trajectory 
is the same in position and momentum and is minimal (that is to say the 
Heisenberg inequality is saturated). Therefore, as announced, the coherent 
state is indeed the "most classical" state. It is also interesting to give the 
explicit form of the wave- function. It reads 

/ , \ 1/4 

f _ 1 c ,-ikt/2 iqp cl (t)-k[q-q cl (t)] 2 /2 



& a (q,t) = I -J e -*'W-' l '-* 1 W | l\ (174) 

where the phase factor is defined by e* e ° = e^ Q ~ a 

The above expression is defined in real space. However, if one wants to 
follow the evolution of the system in phase space, it is interesting to introduce 
the Wigner function [T51 H5J [SSI EZ] defined by the expression (for a one- 
dimensional system) 

1 f . -u. / u \ . , _ / u 



W(q,p,t) = — duV*(q- e-^(g+ \,t) . (175) 
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Fig. 4. Wigner function (|176[) for the coherent state |a) at different times. The (ar- 
bitrary) values qo — 1, Po — 1 and k — 2 have been used for this figure. This implies 
a = V2e i7r/4 and, see Eqs. ([165]) . p d = 2^2 sin (tt/4 - t) and q cl = ^2 cos (tt/4 - 
The upper left panel represents the Wigner function (|176[) at time t = while the up- 
per right, lower right and lower left panels correspond to W(p, q, t) at time t = tt/2, 
t = tv and i = 37r/2 respectively. The wave packet follows the periodic (ellipsoidal) 
classical trajectory in phase space and its shape remains unchanged during the mo- 
tion. 

A system behaves classically if the Wigner function is positive-definite since, 
in this case, it can be interpreted as a classical distribution. In addition, if 
the Wigner function is localized in phase space over a small region corre- 
sponding to the classical position and momentum, then the corresponding 
quantum predictions become indistinguishable from their classical counter- 
parts and we can indeed state that the system has "classicalized" . For the 
wave- function (|174[) . the Wigner function can be expressed as 

W(q,p,t) = le-M'-toWlVife-^MJ' . (176) 

7T 

It is represented in Fig. |4j We notice that the Wigner function is always 
positive and, therefore, according to the above considerations, the system can 
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be considered as classical. Moreover, W(p, q, t) is peaked over a small region 
in phase space and the wave packet follows exactly the classical trajectory (an 
ellipse), as is also clear from Eq. (|176|) . This confirms our interpretation of 
the coherent state as the most classical state. 

To conclude this sub-Section, let us notice that the coherent state \a) can 
be obtained by applying the following unitary operator on the vacuum state 

| a) =c Qat - Q * a |0). (177) 

This equation should be compared to Eqs. lj90|) and (|9l"j) . We see that the 
argument of the exponential is linear in the creation and annihilation operators 
while it was quadratic in the case of the squeezed state. 



5.2 Wigner Function of the Cosmological Perturbations 



In order to study whether the (super- Hubble) cosmological perturbations have 
"classicalized" , we now use the technical tool of the Wigner function intro- 
duced before. The first application to cosmological perturbations was made 
in Refs. [12, 33 . For a two-dimensional system (here, we have in mind fi^ and 
fj, k for a fixed mode fc), the generalization of Eq. ()1 75|1 is straightforward and 
reads 



* //d.d^*(^-f,Mi-i 



(2 

/ u v 



-ipgu—ip'.V 



(178) 



where the wave-function is given by the expressions (I149| . Since we have 
to deal with Gaussian integrations only, the above Wigner function can be 
calculated exactly. One obtains 



1 



2n$tf2 k 
1 



exp 



x exp 



(^+23/4/4) 



(179) 



It is represented in Fig. \5\ The first remark is that the Wigner function is 
positive (as expected since we deal with a Gaussian state) and, therefore, can 
be interpreted as a classical distribution. However, as shown in Fig. [5l and 
contrary to the case of a coherent state, W is not peaked over a small region 
of phase space. We are interested in the behavior of the Wigner function for 
modes of astrophysical interest today. These modes have spent time outside 
the Hubble radius during inflation and, as a consequence, their squeezing 
parameter r is big. Therefore, it is convenient to express f2 k in terms of the 
squeezing parameters 
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Fig. 5. Wigner function of cosmological perturbations obtained from Eq. (|179[) (for 
a one-dimensional system). The squeezing parameter r is chosen to be r — 0.1, 
r = 0.5, r = 1 and r — 2 for the left upper, right upper, left lower and right lower 
panels respectively (it is not the same time ordering as in Fig. [4] because, in the 
present case, the motion is not periodic). The other squeezing parameters are taken 
to be <j> — 7r/6 and 8 — 0. As can be noticed in this figure, the Wigner function 
remains positive. Since the squeezing parameter increases with time, the different 
panels correspond in fact to the Wigner function at different times. At initial time, 
the quantum state is the vacuum and, therefore, the Wigner function is that of a 
coherent state, compare the left upper panel with Fig. [3] Then, the Wigner function 
develops the "Dirac function behavior" discussed in the text that clearly appears on 
this plot. 



k cosh r — e sinh r 
k 2 coshr + e -2 '^ sinhr ' 



and to take the strong squeezing limit, r — *■ +oo. One has 



Rrtfc-O, 3^ fc ^^. (181) 

2 cos ffl 



This implies 
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where (5 denotes the Dirac function. The above limit is clearly visible in Fig. [5] 
for r > 1 (lower panels). 

Therefore, in the large squeezing limit, the Wigner function is elongated 
along a very thin ellipse in phase space. At first sight, this means that the 
system is not classical since one can not single out a small cell around some 
classical values (/i c i,Pci) that would follow a classical trajectory as it was dis- 
cussed before. On the other hand, as already mentioned above, the Wigner 
function remains positive. This means that the interference term which makes 
the system quantum in the sense that the amplitudes rather than the proba- 
bilities should be summed up have become negligible. Therefore, in this sense, 
the system is classical or, more precisely, is in fact equivalent to a classical 
stochastic process with a Gaussian distribution (given by the term ^<l r *). We 
see that the nature of this classical limit is quite different to what happens in 
the case of a coherent state: we cannot predict a definite correlation between 
position and momentum but we can describe the system in terms of a classical 
random variable. In practice, this is what is done by astrophysicists: in par- 
ticular, the quantity at m in Eq. (|151|) is always treated as a Gaussian random 
variable and any detailed quantum-mechanical considerations avoided. 

As argued in Ref. [13], the system has become classical (in the sense ex- 
plained before) without any need to take into account its interaction with the 
environment. This is "decoherence without decoherence" as stressed in the 
above-refered article. More on this subtle issue can be found in Refs. [13], 138] . 
Of course, the question of whether the wave-function of the perturbation has 
collapsed or not (and the question of whether this question is meaningful in 
the present context and/or dependent on the interpretation of quantum me- 
chanics that one chooses to consider) is even more delicate |63j and we will 
not touch upon this issue here. 



5.3 Wigner function of a Free Particle 

We now consider the case of a free particle since it shares common points 
with the case of cosmological perturbations as first noticed in Ref. [38]. At 
the beginning of this section, we mentioned that, in the particular case of 
the harmonic oscillator V(q) oc g 2 , the quantum mean value of the position 
and momentum operators always follow the classical trajectory whatever the 
quantum state |^) in which the system is placed. There is obviously another 
situation where this is also the case: the free massive particle where V(q) = 0. 
The wave- function is given by (we take m — 1) 

1PY„ +\ — f 2(l2 \ 1 1 -it a .n- 1 (2t/a 2 )/2+ik (q-q )-iklt/2 

xe -(<?-<?0-M) 2 / (a 2 +2»t) ^ (183) 
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Fig. 6. Wigner function (|188|) of a free particle. The parameters chosen are feo = 1, 
go = 1 and a — 1 which is in fact equivalent to considering the dimensionless 
quantities ap, q/a and t/a 2 . The left upper panel corresponds to t/a 2 — while the 
right upper, left lower and right lower panels represent the Wigner function at times 
t/a 2 = 0.8, 2, 3 respectively (here, the time ordering is similar to that in Fig. [S]and, 
therefore, different from that in Fig. [4}- Initially, the wave packet is well-localized 
in phase space and, as time goes on, the spreading of the wave packet becomes 
apparent. 

where we have assumed to the wave packet is centered at q — qo at t — to . The 
parameter a represents the width of the wave packet while fco parameterizes 
its velocity. The means of p and q can be expressed as 

(p) = fc = Pd , (q) = qa + k t = q d . (184) 

Therefore, as announced, they follow exactly the classical trajectory. But, as 
argued in the sub-Section devoted to coherent states, one must also compute 
the dispersions. Straightforward calculations lead to 

(f) = k 2 Q + 1 , (a 2 ) = (a 4 + At 2 ) + (q + k tf , (185) 
from which one deduces that 
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Ap=~, Aq=^Jl+4^. (186) 
a I V a* 

At t = 0, one has AqAp =1/2 and the Heisenberg bound is saturated: at 
initial time, the wave packet is minimal. But then, and contrary to the case 
of the harmonic oscillator, the dispersion on the position is increasing with 
time (while the dispersion on the momentum remains constant). The wave 
packet does not keep its shape unchanged while moving as it was the case for 
a potential V(q) cx q 2 . This is the well-known phenomenon dubbed "spreading 
of the wave packet" . However, it exists a quasi-classical interpretation of this 
situation. Indeed, when i>o 2 , one has 

Aq ~ - = Apt = Av c it , (187) 
a 

which reproduces the classical motion. On the contrary, for small times, Aq 
must take values very different from the classical ones in order to satisfy the 
Heisenberg inequality. Therefore, in the regime t ^> a 2 , the system is classical 
but in a sense slightly different from the one encountered in the harmonic 
oscillator case. 

Let us now calculate the Wigner function of the free particle. One obtains 



W{p, q,t) = - exp 

7T 



x exp 



2a 2 



a 



4 +4t 2 



(q - q c \f 



a 4 + 4t 2 r ■ 1 



2a 2 



p - pci -^f^ {q - qcl) 



(188) 



This equation is similar to Eq. (|179p . The Wigner function is the product of 
one exponential factor whose argument is proportional to " [a — q 2 ^" [hid- 
den in the term il/*^ in Eq. (|179[) ] and of another exponential term whose 
argument has the form "[p — p c i — f(t) (q — q c i)] , where f(t) is a function of 
time only. Therefore, the classical limit of cosmological perturbations can also 
be understood in terms of the (quasi-) classical limit of a free particle, as 
discussed in the previous paragraph. 

The Wigner function (|188| is represented in Fig. [6l This plot confirms 
the interpretation presented above. First of all, the Wigner function remains 
positive which indicates that a classical interpretation is meaningful. At initial 
time, the Wigner function is well-localized because the wave packet is minimal. 
Then, at time goes on, the spreading of the wave packet causes the spreading of 
the Wigner function in phase space. Clearly, this case bears some ressemblence 
with that of cosmological perturbations, compare Figs. [5] and [5] Therefore, 
inflationary fluctuations on large scales become classical in the same sense 
that a free particle is classical far from the origin. A much more detailed 
description of this analogy can be found in Ref. [38] , 
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6 Conclusions 

In this review, we have presented a pedagogical introduction to the theory of 
inflationary cosmological perturbations of quantum-mechanical origin, focus- 
ing mainly on its fundamental aspects. We have shown that the mechanism 
responsible for the production of the initial fluctuations in the early universe 
is in fact similar to a well-known effect in quantum field theory, namely the 
Schwinger effect. It is indeed the "interaction" of the quantum perturbed met- 
ric 5g^ v (and of the perturbed inflaton field Sip) with the background gravi- 
tational field which is responsible for the amplification of the initial vacuum 
fluctuations in the same way that pair creation can occur in an external elec- 
trical field. Because the gravitational field in the early universe can be strong 
(in Planck units), this mechanism leads to observable effects, in particular to 
CMB temperature fluctuations. Therefore, these fluctuations originate from 
a remarkable interplay between general relativity and quantum mechanics. 

There is also another aspect associated with the inflationary mechanism 
discussed above that could be relevant to probe fundamental physics. In a 
typical model of inflation, the total number of e-folds is such that the scales 
of astrophysical interest today were initially not only smaller than the Hubble 
radius but also smaller than the Planck (or string) scale. As a consequence, 
the WKB initial conditions discussed before are maybe modified by quantum 
gravity (stringy) effects and, therefore, this opens up the possibility to probe 
these effects through observations of the CMB [68l [69] . Of course, an open 
issue is the fact that it is difficult to calculate how quantum gravity/string 
theory will affect the initial conditions. Nevertheless, it is possible to draw 
some generic conclusions. Firstly, the standard initial condition consists in 
choosing only one WKB branch. Hence, any modification amounts to con- 
sidering that the second branch is present. As a consequence, super-imposed 
oscillations in the power spectrum unavoidably appear, the amplitude and the 
frequency of these oscillations being unfortunately model dependent [68j (for 
the observational status of these oscillations, see Refs. [69],I61|). Secondly, the 
presence of the second WKB branch means, in some sense, the presence of 
particles in the initial state and, therefore, there is potentially a back-reaction 
problem. Generically, the larger the amplitude of the super-imposed oscilla- 
tions is, the more severe the back-reaction issue. On the other hand, predicting 
the effect of the energy density of the perturbations is difficult and it is not 
clear whether this will spoil inflation or, for instance, just renormalize the 
cosmological constant [701 HI] ■ These problems are still open questions but it 
is interesting to notice that the inflationary scenario is rich enough to provide 
yet another means to learn about fundamental physics. 
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